A Rediatively Light Stop Saves the Best Global Fit for Higgs Boson Mass and Decays 
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The LHC discovered the Standard Model (SM) like Higgs boson with mass around 125 GeV. 
However, there exist hints of deviations from Higgs decays. Including the Tevatron data, the devi- 
ations can be explained by the extremely mixed stop sector in the sense of best global fit (BGF). 
We analyze the relations among the competing reduced coupling hGG, Higgs boson mass, and LHC 
stop mass mj- lower bound at the tree- and one-loop level. In particular, we point out that we use 
the light stop running mass in the Higgs boson mass calculation while the light stop pole mass in 
the Higgs decays. So the gluino radiative correction on the light stop mass plays the crucial role. 
Its large negative correction saves the GBF in the Minimal Supersymmetric SM (MSSM) and the 
next to the MSSM (NMSSM) constrained by the perturbativity. Moreover, a light stop is predicted: 



in the MSSM if we set the gluino mass M3 < 4 TeV, we have 
imposing the least tuning condition we get mr ~ 130 GeV for M3 : 



< m t ; while in the NMSSM, 
1.5 TeV. 
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Introduction: The LHC has not found any direct signa- 
ture of supersymmetry (SUSY) so far. But at the 5a level 
the discovery of the Standard Model (SM)-likc Higgs bo- 
son with mass mh —125 GeV [1, 2], indirectly supports 
SUSY from several aspects. In the first, if this spin-0 
particle is confirmed as a fundamental particle, we must 
explain the corresponding quardratic divergency prob- 
lem, and SUSY is the most elegant solution. Next, one 
may wonder why the Higgs boson is so light, i.e., why is 
it not 150 GeV or even heavier? But this is not a puzzle 
in the SUSY models, which actually are suffering from 
the opposite puzzle: why is the Higgs boson so heavy? 

Last but never the least, the discovered Higgs boson 
is not exactly SM-like and shows some notable devia- 
tions which may imply new physics. Both the gluon 
fusion (GF) and vector boson fusion (VBF) produc- 
tion modes show significant di-photon excess. For the 
pp — > h — > V*V channel there exists a mild deficit in 
the CMS experiment [2] but small excess in the ATLAS 
experiment [1]. Including the Tevatron data, the best 
global fit (BGF) gives [3] 
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with all the other couplings fixed. The reduced couplings 
c gn are sensitive to colored and/or charged particles be- 
yond the SM, and thus indicating new physics. 

In SUSY models, stau can increase c 7 but does not 
affect c g [6] while the stops, which carry the SU(3)c 
and U(1)em quantum numbers simultaneously, may al- 
ter the reduced couplings towards Eq. (1). Practically, 
if the stop-loop interferes destructively with the top-loop 
contributing to h — > gg while constructively with the 
W— loop contributing to h —> 77, then the stops can im- 
prove the global fit by an amount A% 2 ~ 10 (compared 
to the SM fit). However, the above BGF in Eq. (1) can 
not be realized generically in the SUSY models due to 
the relatively large SM-like Higgs boson mass. In this 



letter, we investigate the possibility of accommodating 
such BGF in the conventional minimal supersymmetric 
SM (MSSM) and the next-to-the MSSM (NMSSM) [4]. 
We find that the BGF can be obtained if the gluino ra- 
diative correction to the light stop pole mass is negative 
and large. 

c g competing with rrih and mj : The stop sector is the 
focus of SUSY phenomenology right now due to various 
reasons. Firstly, the stop corrections [7] lift mj above 
the LEP bound and even to 125 GeV indicated by the 
LHC. At the one- loop level, the corrections can be ap- 
proximated as 



(2) 




4ir 2 v 2 



with v = 174.6 GeV, 7] t = X t /mj with X t = A t - H cot (3 
the stop mixing parameter, and mj = ^/to^ the ge- 
ometric mean of two stop masses. Here, A t is the top 
quark trilinear soft term and fi is the Higgs bilinear mass 
term in supcrpotcntial. When rj t <C 1, the mixing term is 
ignorablc and the leading logarithm dominates the cor- 
rections. While in the maximal mixing scenario with 
r\t m V6 [8], the mixing term effect maximizes. However, 
as r\t > 2\/3 this effect begins to act reversely. We will 
find that it is the typical case to approach the BGF. 

Next, the stop sector can change the reduced couplings 
of hGG and /177 simultaneously. In terms of the BGF 
in Eq. (1), the stop-loop should destructively interfere 
with the top-loop and thus modify the total reduced hGG 
coupling to be [9] 
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-0.84, (3) 
where the heavier stop contribution is ignored. Accord- 
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ingly, the reduced /177 coupling is modified to be [9] 



1.28- 0.28c o 
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which is valid for ~ 125 GeV and cy ~ 1. 

Finally, the stops carry colour charges and have a close 
relation with the pattern of electroweak symmetry break- 
ing. Thus, they are among the main search particles 
at the LHC. Presently, the lighter stop is bounded by 
mj- > 130 GeV (details is postponed). Then we get 
mj- > 600 GeV from the approximation which is valid 
for highly mixed stops: 



\X t \^{2C G fl 2 -^m Tlj , (5) 



with Cg ~ 4:(l — c g ) + mf /m~ . We have taken the simpli- 



fying assumption vnq R 



, but it can be abandoned. 



Synthesizing these aspects, we obtain the competing 
relations among c g and mh, mj , which make it difficult 
to realize the BGF in the (N)MSSM. (A) At least in the 
MSSM, to get mh — 125 GeV a large mj and properly 
mixed stop sector are necessary. It means that generi- 
cally the mixing term in Eq. (3) is ignorablc, while the 
lighter stop enhances c g instead of the other way around. 
(B) The NMSSM, in spite of its advantage in lifting mh 
which weakens the relation with the stop sector, is also 
impossible. This is blamed to the lower bound on mj , 
which leads to a big suppressing m\ /m~ -C 1 in c g . Con- 
sequently, if the mixing term can still make c g ~ —0.8, we 
likely need rj t > 12 to compensate for that. In turn, such 
a large mixing term, which grows quarticly in At, will 
decrease m/j. However, this behavior is only an approx- 
imation, and it underestimates Ah in the large mixing 
region. In the numerical plot we will adopt the complete 
expression. 

A radiatively light stop mass as a savior: The previ- 
ous robust relations from the two-parameters stop system 
can be relaxed, if we take into account the sizeable one- 
loop radiative correction to the lighter stop (the much 
smaller correction to the heavier stop is neglected here). 
Such correction is crucial, since the stop mass used to cal- 
culate At is the running DR mass, while the stop mass 
in Eq. (3), used to calculate c g , is the pole mass. They 
are related by the equation 
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where the running mass mj(Q) is defined at the renor- 
malization scale Q = mj . H(p 2 ,Q) is the self-energy, 
dominated by the QCD and top- Yukawa radiative cor- 
rections from the gluon, gluino, stop and chargino loop, 
etc. When the positive self-energy is sufficiently large, 
the pole mass can be much lighter than the running mass. 

The gluino mass plays the primary role to reduce the 
running mass. There exists some uncertainties from the 
Yukawa and two-loop corrections. In the analytical anal- 
ysis we will not go through the details of this discussion, 
but instead introduce an extra parameter oj t . the ratio of 



pole mass and running mass, to account for those. Illus- 
tratively, at the order 0(a s ), u>t is roughly given by [10] 
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The approximation is valid for M^,m^ 2 {Q) 3> mj (Q) ~ 
riif. Note that we have mj 2 > Q. Given a heavy gluino 
(required by the LHC bound), the running mass can read- 
ily be reduced by a factor 2 or 3. As a consequence, one 
can work in the maximal mixing scenario r] 2 (m^ L ) ~ 6, 
while rj 2 (pole) is several times larger than rj 2 {fn^ L ) . In 
this way, the radiatively reduced stop mass offers a solu- 
tion to reconcile the serious tension between mh and c g . 
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FIG. 1: Schematic plots of m h = 125 GeV, c 9 = -0.8 
(red thin solid lines) and mj- = 130 GeV (dashed lines) 
on the m^ L — Xt plane. For comparison, we choose ujt = 
0.5,0.4,0.3,0.2 for c g and m^, respectively from the left to 
right. In plotting mj, we have fixed tan /3 = 25 in the MSSM 
and tan /3 = 2, A = 0.6 in the NMSSM. The running top mass 
is taken to be 155 GeV. 

The schematic analysis in Fig. 1 manifests the point. In 



the MSSM, as w t decreases the c„ 



-0.8 curves, which 



are denoted by the red solid lines, move towards the left 
and meet with the mh = 125 GeV curve at a large mj . 
As for the NMSSM, such meeting still needs w t < 0.6 (we 
will turn back to this point later). In this case, the lower 
bound on the stop pole mass m^ =130 GeV means that 
the running mass is above 220 GeV. In turn, we require 
a rather heavy gluino mass to account for the big gap. 
Additionally, from Fig. 1 one can see that for a given u) t , 
the c g = —0.8 curve intersects with the =130 GeV 
curve at large mj (understood similarly to Eq. (5)). And 
its value increases as Wt decreases. 

Wt < 0.6 is owing to the NMSSM perturbativity up 
to Mqut- It renders A < 0.7 and imposes the tree-level 



3 



bound mh < Avsin/3 ~ HOGeV. Then a moderate stop 
correction and/or pushing effect is still indispensable. In 
turn, in light of the relations between competing mq , rrih 
and c g , the NMSSM with w t ~ 1 fails. But in ASUSY [11] 
with A ~ 2, or other non-decoupling F/D— term mod- 
els [12], tree-level mh can be sufficiently heavy such that 
large rjt{> 12) is endured. 
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FIG. 2: Top figure: The plot of the parameter space for 
mj- -|At|-M3, where we have required 0.7 < c? g < 0.8 and 
123 GeV< m h <127 GeV; Bottom figure: The plot of distri- 
bution of Wt. 

Numerical search in the (N)MSSM: For a more 
precisely quantitative study, we turn to the NMSSM- 
Tools 3.2.0 [13]. Firstly we get the MSSM limit from 
the NMSSM, e.g., setting it to the decoupling limit: 
A = k = 0.01, tan/3 = 15, fj, = 500 GeV, A K = 800 GeV, 
Ma = 1.5 TcV. One may vary them freely, but it will not 
bring much difference to the configuration of the three- 
parameter-space m-j: -\At\-M3, projected on the mj -\At\ 
plane. From Fig. 2, it is seen that the right edge of the 
triangle can extend when we increase M3, etc. Accord- 
ingly, the right edge of the wt distribution slowly moves 
towards the right, where wt becomes smaller. In sum- 
mary, a multi-TeV gluino, \At\ and quite small ut are 
the main features of the parameter space in the MSSM. 

Now we turn our attention to the limit which manifests 
the NMSSM specific effects enhancing m,h [14, 15]. It is 
characterized by A ~ 1, tan/3 ~ 1 and/or a small /z. Es- 
pecially, a small /i, which is consistent with naturalness, 
leads to not only the Higgs doublet-singlet (HS) mixing 
pushing upward mh [14] but also the di-photon excess 
by suppressing Br(/i — > bb) [14, 16]. But in view of the 
BGF, the stop-loop instead of the mixing effect should 



account for the excess. Therefore, we require a small HS 
mixing angle such that hff and hVV are not modified 
much, as forces us to accept some degree of tuning, from 
either a relatively large (i or a mild cancellation in the 
Higgs mass matrix. 

We restrict the discussion to a smaller [i and explore 
how natural the NMSSM can be after requiring the BGF. 
Roughly speaking, this amounts to find the gluino lower 
bound [9, 14], which has been shown to be large. \A t \ 
also should be sufficiently small [9, 14, 17]. We give a 
benchmark point here. The inputs at the SUSY scale 
are M Y = HO GeV, M 2 = 740 GeV, M 3 = 1500 GeV, 
A t = -1530 GeV, and mj L = 588 GeV. Also, tan/3, fi 
and the NMSSM specified parameters are tan/3 = 2.17, 
A = 0.666, k = 0.253, A x = 510 GeV, A K = -248 GeV, 
and \i = 291 GeV. The resulting particle spectra contains 
the SM-likc Higgs, light stop and lightest neutralino with 
the dominant bino component. Their pole masses are 
given by m h = 124.8 GeV, M x o = 100.0 GeV, and = 
130.3 GeV. In particular, c g = —0.84 is generated. 

We would like to make some comments. (A) The stop- 
chargino loop generates flavor transition b — > sj. How- 
ever, its rate is safely small if we take tan/3 < 20 and fi 
a few hundred GeV [18]. (B) Due to the large A t — term, 
the split stop-sbottom may render Ap orders larger than 
the experimentally allowed bound < 10~ 3 . Small mj- is 
favored to lower the size of Ap, and some cancellation 
from is needed. (C) We can relax the lower bound 
on sfightly, e.g., to 125 GeV, by means of introduc- 
ing R— parity violation operators U C D C D C . It makes the 
very light stop missed at the LHC due to the absence of 
missing transverse energy. In this case, we do not need 
the heavy gluino and the BGF can be naturally realized 
in the NMSSM. 

A light stop escapes from the LHC: In both models, 
if M 3 < 4 TeV, we predict < m t , see Fig. 2. As 
long as wc impose the minimal fine-tuning condition, it 
is found that the NMSSM allows M 3 < 1.5 TcV and gives 
m+ ~ 130 GeV which will be checked safe. 

The LHC search for the stop is sensitive to its decay 
topology and the particle spectrum. As a case in point, 
we consider the benchmark point in the NMSSM. It con- 
tains the relevant particles, i.e., the light stop with mass 
130 GeV and the 100 GeV bino-likc neutralino LSP Xv 
Due to the kinematics under consideration, t\ dominantly 
(nearly 100%) decays to c and Xi y i a the- loop-induced 
process. The process at the LHC is 



99 -> tit*i 



c + c + 2 X ?. 



(8) 



In this case, the monojet and jets+i?™ lss signatures give 
the strongest bounds [19]. Actually, the monojet bound 



is effective only when < 



< 0.2?™, 



Oth- 



erwise, a putative case in this work, it is ineffective and 
need the inspection to jets+i?™ lss , discussed in the fol- 
lowing. 

At first, we use MadGraph5 vl.4.7 [20] to generate the 
stop pair, accompanied by (< 3) additional initial state 
radiative (ISR) jets. In calculating this cross section, we 
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« + jj 


tt + (> 3j) 


Cross section (pb) 


42.08 


15.73 


3.84 


0.87 



take the modified MLM matching with xqcut = 30 GeV 
and showcr-kt scheme with QCUT=50 GeV. Then after 
matching the LO production cross sections of t\t\ + nj 
at the 7 TeV LHC, are displayed in the following: 



(9) 



Moreover, the k-factor calculated by the Prospino2 [21] is 
about 1.5 for mj- ~ 130 GeV. After these, we use Pythia 
to do the particle decay, parton showering and hadroniza- 
tion [22]. PGS is chosen to simulate the detector effect. 
Finally, employing the cuts adopted in Rcf [23], we find 
that only the events with two or more ISR jets survive. 
The most sensitive signal regions and the corresponding 
event numbers (normalized to 4.7/fb) are given by: 



Signal Region 



Number of Event 



SRC-med 



2.03(18) 



SRC-loo 



15.10(58) 



SRE-med 



1.02(12) 



SRE-loo 



4.07(84) 



where the numbers in the brackets are the observed upper 
limits on the event numbers beyond the SM. Clearly, the 
present ATLAS jets+i?™ lss search can not exclude our 
light stop scenario. We are awaiting the ongoing LHC to 
probe such scenarios. 

Conclusions and discussions: Stops are potential to 
accommodate the BGF. At first glances, due to the re- 
lations among the competing c g , mj- , and m^, there is 



merely a sporting chance to realize the BGF in the MSSM 
and NMSSM. Fortunately, the gluino radiative correc- 
tion on the light stop mass can make its pole mass much 
lighter than the running mass, and then the BGF can be 
realized. If we reasonably let M3 and \A t \ below 4 TeV, 
a light stop typically below the top mass is predicted. 
In particular, in the NMSSM we explore the less fine- 
tuned region to the most extent, and find that the model 
allows M 3 ~ 1.5 TeV as well as m 7i tn 130 GeV. Such 
light stop is not excluded by the ATLAS latest data since 
the decay t\ -> c + Xi with x° about 100 GeV. We en- 
courage the LHC experiments to tailor a search for the 
light stop inspired by the Higgs BGF. Moreover, besides 
the light stop, the lightest neutralino, which may be an 
successful dark matter candidate, also deserves further 
attention. In our scenario, its mass falls into the region 
ttlt > M v o > 0.83mr , i.e., between 100 GeV and the 

ti Xi ti> ' 

m t . We leave it for a future study. Finally, after the com- 
pletion of this work we are aware of [24] , which gives a 
stringent constraint on the BGF inspired stop sector from 
the vacuum stabilities. But their constraint is based on a 
rather simplified analysis and only valid at the tree-level. 
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